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We review recent work in the lattice approach to random surfaces and quantum gravity. Our task is made 
somewhat easier by some very interesting results, particularly in four dimensions, that have appeared recently and 
which are reported elsewhere in these proceedings. Inevitably, given the scope of the review and the limitations 
of space, the presentation will omit work of importance and be telegraphic in discussing work that is included, 
for which apologies are offered in advance. After the customary brief historical introduction we work our way 
in dimensional order from one up to four dimensions before closing with some remarks on the relation, if any, 
between the various lattice models and "real" 4D gravity. 



1. INTRODUCTION 

It is probably fair to say that lattice gauge the- 
ory in the large is now a mature subject. The 
theoretical underpinnings are understood and the 
numerics are under control. The subject of this 
review is a much more speculative application of 
lattice methods, namely to the study of quantum 
gravity in various dimensions. Geometry itself be- 
comes dynamical in gravitation, so it is clear that 
unlike standard lattice theories we will be dealing 
with dynamical lattices of some sort in which the 
geometry and the matter living on the lattice, if 
any, are in interaction. 

The first question to ask before embarking on 
the numerical investigation of any theory is why 
bother? In the case of gravity, the answer is quite 
clear: the Einstein-Hilbert action for general rel- 
ativity is perturbatively non-renormalizable, es- 
sentially because of the dimensionful coupling. 
Although this does not exclude a valid theory |^, 
it does force us to employ methods other than 
those of perturbative quantum field theory to in- 
vestigate the model. There are further problems 
lurking even after the rotation to Euclidean sig- 
nature which is generically necessary to obtain 
tractable simulations on the lattice. These are 
due to the unboundedness below of the Euclidean 



'^An example of a perturbatively non-renormalizable tlie- 
ory that does have a well-defined continuum limit is the 
4D 0(4) non-linear sigma model 



action, due to conformal mode fluctuations, 
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where M is our spacetime manifold. At first sight 
this might seem to render even the Euclidean par- 
tition function 
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Dgab] exp {-S) , 



ill-defined. We might be saved, however, by the 
measure in the path integral giving negligible 
weights to such configurations. In order to see if 
such a happy occurence does take place we need 
the framework of some regulated lattice theory 
where the integral/sum over configurations can 
be clearly defined. 

It is perhaps worth strengthening the case for 
numerical investigations by noting that minimal 
tinkerings with the quantum field theory don't 
work. Introducing higher derivative terms gives 
a renormalizable theory with action 
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which is is even asymptotically free for suitably 
chosen a and f3 Q but, as might be expected, 
one runs into problems with unitarity. More rad- 
ical solutions, SUSY, SUGRA, strings, p-branes, 
M-theory, non-commutative differential geome- 
try have all generated a lot of beautiful math- 
ematics (and even more preprints) but getting 
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back to boring old four dimensions and a non- 
supersymmetric, or at least broken supersymmet- 
ric particle spectrum, poses problems of varying 
degrees for all of them. 

A more modest line of attack is to attempt 
to formulate a non-perturbative approach to the 
Einstein-Hilbert action and its close relatives. 
There are two complementary schemes that have 
been used in most recent work, Regge Calculus 
(RC) and Dynamical Triangulations (DT), which 
we now describe briefly in turn. 

Regge calculus (RC): was originally sug- 
gested by Regge in 1961 |^ as a means of dis- 
cretizing classical GR. The recipe for discretizing 
the D-dimensional spacetime manifold M was to 
replace it by D-simplices. Curvature then resides 
on the D — 2 dimensional hinges between the sim- 
plices (i.e. vertices in 2D gravity, links in 3D 
gravity and faces in 4D gravity). In the RC ap- 
proach the discretized equivalent of the integra- 
tion over metrics is usually implemented by tak- 
ing a fixed connectivity simplicial complex with 
varying edge lengths. Regge and Ponzano came 
back for a second bite with a theory of quantum 
3D gravity in 1968 and more recently there 
has been much numerical work on the quantum 
theory in AD by Hamber Q and the Vienna group 

Dynamical Triangulations (DT): The re- 
cent work with this approach is an outgrowth of 
ideas that had their genesis in the theory of ran- 
dom surfaces and 2D gravity Q but this variant 
also has a long pedigree, in that the basic ideas 
were formulated Weingarten in 1982 using hy- 
percubes. As with RC, one replaces M by D- 
simplices, but now we take varying connectivity 
in our simplicial complex and fixed edge lengths 
in order to implement a discretized version of the 
integral over metrics. In two dimensions this for- 
mulation is essentially a recipe for the numer- 
ical evaluation of matrix model partition func- 
tions and has been an unqualified success. As a 
cautionary remark it should be noted that one 
does not have the same theoretical underpinning 
in three and four dimensions. 

From the statistical mechanical viewpoint in 
both RC and DT we are dealing with theories 
with a particular sort of annealed disorder com- 



ing from the geometric fluctuations. This is as we 
might expect - the dynamics of the lattice should 
be taking place on the same timescale as any mat- 
ter, so we have annealed rather than quenched av- 
erages to contend with. In RC we have annealed 
bond length disorder whereas in DT we have an- 
nealed connectivity disorder. 

In both RC and DT the discretized actions are 
appealingly simple, being for Regge Calculus in 
4D 

t s 

where At is the area of triangular faces t, Vs is 
the volume of simplices s and the deficit angle is 
given by 

<5t = 27r - ^ es,t 

tes 

and for dynamical triangulations in various di- 
mensions 

S = KiNi-K2N2 (41?) 

S = K37V3-«iA^i (3Z?) 
S = K2N2 {2D) 

where the Nk are the numbers of fc-simplices and 
the various k's are the coupling constants. 

There are various choices to be made in the two 
schemes. In RC the measure for the integration 
over the edge lengths must be chosen according 
to some prescription and a fixed triangulation de- 
cided upon. In DT one must decide which class 
of simplicial complexes to employ - for instance 
in 2D the finite size effects are minimized, rather 
counterintuitively, by allowing various degenerate 
triangulations. In both schemes one must decide 
which topology to fix or whether, indeed, to let it 
vary. 

In DT the moves that change the connectivity 
in D dimensions can be obtained from "slicing" 

+ 1 dimensional simplices appropriately. In 
two dimensions there is a "fiip" move that pre- 
serves the number of vertices in a triangulation 
as well as an insert/delete move, so it is possible 
to cary out either canonical fixed number of nodes 
simulations or grand canonical varying number of 
nodes simulations. The moves in two dimensions 
are shown overleaf: 
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This option is no longer open in three and four 
dimensions where the moves change the number 
of simpHces. These moves are represented in their 
simplest form as so-called (fc, I) moves ||^ which 
take k subsimplices to I subsimplices, and may 
be shown to be ergodic. The moves in the three 
dimensional case are shown below: 




being (1,4) and (2,3) moves respectively. My 
very limited skills as an illustrator defeat me in 
the four dimensional case where the (fc, I) moves 
are (1,5), (2,4) and (3,3). In addition cluster- 
like moves have also been defined, which have the 
effect of cutting off large chunks of the simplicial 
complex and gluing them back on elsewhere. 

The partition function of interest in the AD 
simulations carried out using the moves described 
above is 

Z{k2,k4)= ^ Z7V2,iV4 exp (-5*) 

where 

T{N2M4) 



with W{T) being the symmetry factor for a sim- 
plicial complex and S being the 4_D DT action. 
To discuss possible critical behaviour it is conve- 
nient to define a fixed volume partition function 

Z(iV4,K2) = exp(K2^2) 
T(Ni) 

SO 

Z(k2, K4) = ^ Z{N4, K2) exp {-K4N4) 

As we have noted it is impossible to carry out a 
fixed volume simulation in AD because an ergodic 
set of (fc, I) moves entails changes in the number of 
simplices, so we perform a fudge to stay "close" to 
given A^4 by modifying the action with a gaussian 
term 

S = K4,N4, - K2N2 + 7(Af4 - Vf 

The observables in our simulations are the dis- 
cretized versions of curvature and the associated 
susceptibihties 

XR ^ {Ni) {<R^>-<R >2) 

as well as other geometrical properties such mean 
geodesic distances, hausdorff dimensions and cur- 
vature correlators. 

A schematic phase diagram for the theory is 
shown below, where moving along the bold line 
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representing the infinite volume limit we may 
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reach a critical point at some value of K2 at which 
correlation lengths may diverge and we can hope 
to take a continuum limit. 

Extensive simulations by various groups both 
in DT and RC |||,^ over the past few years 
suggested the following picture: 

• In DT one sees a low K2 "crumpled" phase 
and a large K2 branched polymer like phase 
in both 3D and AD. 

• A fractal baby universe structure was man- 
ifest in iD and AD (and also in 2D) with 
DT. 

• With DT in 'iD there was a first order tran- 
sition between the crumpled and branched 
polymer phases. 

• With DT in AD the transition was second 
order. 

• RC saw a similar phase structure (but not 
necessarily order of transitions). 

• Issues of an exponential bound on the num- 
ber of configurations [pT| and computability 
in DT appeared to have been settled to (al- 
most!) everyone's satisfaction. 

The schema above for DT has the great merit of 
conceptual simplicity - starting with a very sim- 
ple action in both 3-D and AD one has found no 
continuum limit in 3D (fine, there is no graviton 
in 3D) but a continuum limit in 4D (also fine, 
we want to get a graviton here). It has, unfor- 
tunately, the great demerit of apparently being 
wrong, as the latest numerical work presented in 
these proceedings indicates. 

In the rest of this article we will review the re- 
cent work in the field in one, two, three and four 
dimensions. The lower dimensional work is in- 
cluded not just for completeness, a lot of it has a 
direct bearing on the understanding of the 4D 
theory. The ID work on polymers, of course, 
stands on its own merits as does the 2D work 
on random surfaces which have a wide applica- 
tion outside the rather esoteric world of quantum 
gravity simulations and string theory. 



2. ID STRUCTURES AND POLYMERS 

It seems to be a general principle that models 
of dynamical geometry in higher dimensions will, 
given half a chance, collapse to lower dimensional 
structures - the large K2 branched polymer like 
phase in 3D and 4D DT being a prime example. 
It is therefore of some interest to study simplified 
models of branched polymer like objects in their 
own right. Indeed, it has been suggested that in 
4D DT the entire branched polymer phase may 
be critical because of power law correlations ob- 
served there 




A note of caution on these observations has been 
sounded recently by Bialas who considered 
a model of planar rooted planted trees as shown 
above. The partition function for these is just 

Trees 

where the weight, p{T), of a given tree depends 
on how many vertices Uk of order k are present 

with the tfe's being the weights for order k ver- 
tices. 

He considered correlations of the form 
G{fi,t,r) = ^exp(-^7i)p(T) 

T 

where /i is the "cosmological constant" with a 
marked point a distance r + 1 from the root, more 
specifically Root- Tail correlations 

G(/i, t; r) = ^ d(«i)d(z;,+i)e-^>(r) 

Tr 
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where the d{v) are the degrees of the vertices in 
question. He found, with some subtleties in the 
definitions of correlators, G — 0. However, in 
a canonical (fixed number of vertices) ensemble 
there were negative l/r^ correlations due to 
the absence of short trees. The moral is that the 
almost fixed volume constraint in AD gravity sim- 
ulations might be responsible for the power law 
behaviour seen in the branched polymer phase 
rather than being a sign of true critical behaviour. 
The issue awaits further clarification. 

Planar rooted planted trees were also used by 
Ambj0rn et.al. in an analytical implementation 
of fractal blocking . Various MCRG schemes 
have been proposed for two and higher dimen- 
sional DT simulations. The new feature here, 
compared with MCRG on a regular lattice, is that 
the geometry is both irregular and dynamical so 
the choice of basic blocking move for the lattice 
itself is not immediately apparent. One possibil- 
ity is the so-called fractal blocking or baby uni- 
verse renormalization group which makes use of 
the fractal structure of the DT simplicial mani- 
folds ITslJl^ ]. In all dimensions these contain re- 
gions where the lattice necks down, which can be 
visualized as baby universes budding off a mother 
universe. This can be used to define a block- 
ing move as shown above right - extremal baby 
universes (those with no babies themselves) are 
lopped off and the couplings rescaled. 

The simplicity of the tree model allows an 
analytical implementation of this program and 
shows, surprisingly that fractal blocking moves 
you away from a branched structure directly to 
a linear chain. Whether this signifies a possible 
pathology of the blocking scheme or is merely a 
quirk of the model is not entirely clear. The fact 
that fractal blocking gives good results numeri- 
cally in higher dimensions argues for the latter. 

The last result I want to mention in ID con- 
cerns planar rooted trees 0] where the restric- 
tion on having only one branch emerging from the 
root node is relaxed. This technical change gives 
a soluble model 

Z = ^ cxp {-iin{T)) exp (-/?i?(T)) 

T 

where n{T) is the number of vertices, E{T) = 




X]„lnfc(w) and k{v) is the vertex degree. There 
are two phases: a large (3 branched polymer phase 
with lots of short "bushes" ; and a small /3 elon- 
gated branched polymer phase. The transition 
between these is fourth order 

<^-</.o = (^i-^io)'-^™ 

where the free energy is = — In and ^string = 
0.3237... The result is interesting because it pro- 
vides a counterexample to the belief that only 
values of ^string = '^jn for n = 2, 3,4, . . . were 
possible in branched polymer models. Indeed, 
by tinkering with the vertex weights a range of 
"/string valucs Can bc extracted from the model. 

3. 2D RANDOM SURFACES 

We make an extended stop in 2D also as the DT 
approach has its genesis here. The usual model 
of interest is the Polyakov String 

where is the extrinsic curvature squared and 
acts as a stiffness or rigidity term. This discretizes 
neatly to 

S-\Y.(^'- + ^E(l - "A. • "A,) 

<ij> A 
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where the are normals to triangles. On a dy- 
namical triangulation this action can serve as a 
model for a fluid surface, and on fixed lattice (so 
there is no gravity at all) as a model for a poly- 
merized or crystalline surface. Analytically the 
fluid model appears to be always in a rough phase 
but simulations provide evidence for a crumpling 
transition as A is increased. 

More recently an alternative action also based 
on geometrical ideas, the gonihedric string, has 
been proposed by Savvidy et.al. Q 
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<v> 

where d{aij) — \ti — Q^ijl^j C < 1 ^-iid the are 
the angles between neighbouring triangles. 

Let us start by discussing some recent results 
concerning the MCRG on crystalline surfaces. 
There has been some confusion in the past con- 
cerning the nature of the phase diagram and the 
values of the critical exponents for the crumpling 
transition on crystalline surfaces (which it is gen- 
erally agreed does exist). Espriu et.al. ||l9| have 
developed a Fourier accelerated Langevin equa- 
tion approach to the MCRG on crystalline sur- 
faces 

dS 

4>{p,tn+i) = (t>{p,tn) - Me{p)T- 



d(t){y,tn) 



-^Ate{p)7]{p, tn) 



,{A(A + m2)} 



A(A -I- m^) 

where denotes a Fourier transform. They 
blocked 9 geometric operators and measured the 
coupling constant flow, flnding the exponent rj ~ 
1.527(36) in agreement with direct simulations, 
but with a lot less effort. In addition, qualita- 
tive examination of the /3-function gave support 
to the presence of a crumpling transition. 

The aim of showing conclusively that a flat 
phase does exist in the crystalline model (and 
hence a crumpling transition) was pursued by the 
Syracuse group |20|. The answer is not immedi- 
ately obvious because considering an action of the 
form 



S = 



<i]> 



X,-XA-a 



where the a represents an intrinsic length for each 
edge gives an effective action for the stress-strain 
tensor Uap 



VS 2 , V3 2 2 

—a Uapupa + —a u^^ 



in which the elastic constants apparently vanish 
as a ^ 0, so it is difficult to see how a fiat phase 
might be stabilized. 

However, measurement of the various scaling 
exponents for the running bending stiffness (k in 
the notation of A elsewhere in this review!) 
and the running clastic "constants" , /i and A 

K{q) ^ q-\ fl{q) ^ \{q) ~ g"" 

as well as the roughness exponent 

< >- 

all of which are related by scaling relations C = 
1 — |, 1] = 1 — ^ gave ?7„ = 0.50(1) (so from 
the scahng relation rj = 0.75), C = 0.64(2) (thus 
7] = 0.72 from the other scaling relation). These 
consistent estimates of a positive rj confirm the 
existence of a flat phase. 

Having surreptitiously slipped in a discussion 
of surfaces without gravity, we now return to our 
remit to look at work on dynamically triangulated 
surfaces - in other words, matter coupled to two 
dimensional gravity. One of the most interesting 
recent developments has been the formulation of 
MCRG methods for dynamical lattices. The frac- 
tal blocking scheme has already been described, 
but another scheme, node decimation, has been 
pioneered by the Syracuse group |^l|. The idea 
is again simple: we pick a node, remove it, and 
then retriangulate the gap to make sure that we 
still have a triangulation as shown overleaf. 

In practice it is simpler to pick a node at ran- 
dom, flip the surrounding links until it has only 
three neighbours and then remove it. Work with 
this method has concentrated on 2D simulations, 
looking at the flow of geometric operators such as 



Although their fluctuations diverge. 
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where is the number of neighbours of site i, 
when a measure factor a log qi is varied. The 
Ising model 

S = UiGj 

<ij> 

coupled to 2D gravity has also been investigated 
and measurements of ^string 



Z{N) ~ N'l' 



' exp {^XcN) 



after blocking moves confirm a flow to the Ising 

value ^string = -1/3. 

The other main thrust of recent work in 2D 
gravity has been the clarification of scaling be- 
haviour. This was discussed in some detail by 
Simon Catterall in Lattice95 so I will be 

very brief here. The canonical object of inter- 
est in such measurements for the geometric sector 
of the theory is the distribution of geodesic dis- 
tances, given in a grand canonical ensemble with 
cosmological constant ^ by 

which can be shown to scale as 

G^{R) ~ ly - 1/dH 

where dn is the hausdorff dimension of the tri- 
angulation and F is some scaling function. The 
other scaling exponent 77 is related to "/string by 
"/string = ^^(2 — ?/). Thc number of points in a 



radial shell at distance i? in a microcanonical en- 
semble of volume V can be shown to scale as 



(SiR)) 



V 



R 



and is perhaps a more amenable observable to 
measure. 

Simulations p3[ | have shown that dn = 4 when 
the central charge c of the matter living on the 
triangulation is less than one, so rather surpris- 
ingly the back reaction of the matter does not 
affect dn- As the triangulations degenerate to 
branched polymers for large c it is also known 
that da = 2 in this regime. 

What was unclear until very recently was 
whether a diverging matter correlation length 
could be successfully defined on dynamical tri- 
angulations. Indeed, a toy model of spins on a 
restricted class of triangulations displayed a 
transition but no such divergence. However, sim- 
ulations presented at this conference have 
shown that such a diverging correlation length 
does exist in the Ising and 3-state Potts mod- 
els on the full dynamical triangulations and also 
measured the effect of the matter back reaction 
on the scaling functions - a delicate task as it is 
the nose and tail of the distributions where the 
effects are seen. 

There have also been simulations p6[ of vari- 
ants of the standard action including i?^terms on 
dynamical triangulations with actions of thc form 



(6 



a 



E 

<ij> 



(6 - g.t)(6 - gj) 



A recent analytical solution of dually weighted 
matrix models |^ has shown that i?^ gravity is 
in the same universality class as ordinary gravity, 
which is what is seen in these simulations. This 
is rather similar to particles paths weighted by 
curvature, which behave as random walks unless 
the curvature coupling is taken to infinity. 

One enduring lacuna in our understanding of 
2D gravity is the nature of the so-called c = 1 
barrier. Matrix models and continuum methods 
both fail to get beyond c — 1, but simulations see 
no immediate evidence of pathologies for c val- 
ues not much larger than one. It is only at large 
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c that triangulations collapse to branched poly- 
mers. One simulation that casts some light on 
this problem appeals to methods from electrical 
engineering to probe the nature of the random 
surface . One treats a triangulation of spheri- 
cal topology as a random network of unit resistors 
and puts current in and out at two vertices, mea- 
suring the voltage drop across another two. The 
voltage-current relation is 

V{zi) - V{z2) = R{zi,Z2; Z3, Z4) I 

where the resistance is 



R 



Ztt 



r is the resistivity constant and [zi, Z2, -^s, ^4] is 
the cross-ratio, which appears because of the con- 
formal invariance. If the theory is to have a sensi- 
ble continuum limit the distribution of measured 
r values should become sharper as the lattice size 
is increased, which happens for c < 1 but not, 
apparently, for c > 1. This suggests that the con- 
formal structure of the surface is breaking down 
at c = 1. It is also possible to extract the distri- 
bution of modular parameters on a toroidal trian- 
gulation, but finite size effects are rather stronger 
here, so the evidence for something nasty happen- 
ing just above c = 1 is weaker. 

I have so far relegated RC to a few cursory 
sentences, which I now rectify. Some older simu- 
lations of the Ising model on a 2D RC lattice gave 
an unpleasant surprise ||2g[| , finding results consis- 
tent with the Onsager (flat 2D lattice) exponents 
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rather than the expected KPZ/DDK exponents 
for the Ising model coupled to 2D gravity. If they 
were to be believed RC was failing to simulate 
gravity in 2D. Simulations of the Ising and Potts 
models with DT do find the correct KPZ/DDK 
exponents. 
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An obvious first reaction is that there is a problem 
with the measure for integrating over the edge 
lengths in RC, as this is not prescribed at the 
outset. In general a local measure is used, one 
common choice being 



n 



for various a. It has been pointed out that if one 
thinks of the lattice model as being a discretiza- 
tion of an exact diffeomorphism invariant theory 
then gauge fixing and the resulting Fadeev-Popov 
determinant are obligatory fsot] . The resulting 
non-local action would be very difficult to sim- 
ulate even in 2D and things get worse in higher 
dimensions. An alternative point of view is that 
exact diffeomorphism invariance is recovered in 
the continuum limit even when broken on the lat- 
tice, so we can get away with local measure |3l|] . 

There have been a series of very careful investi- 
gations by Janke and Holm of both pure gravity 
and gravity plus matter in 2D RC with local mea- 
sures [p2| . In the pure gravity case they looked 
at an action 



S ■ 



A-t- 

4 



A 



A 
a 



which has two scaling regimes 

A « 1 

Z{A) - A-^'^*-^ exp (^-^) exp (-AflA - £,A 

giving = 2 - 2(1 - g), Sc = 16^2(i _ ^)2 
the standard 2D gravity regime 

A » 1 

Z{A)^A''^*~^eM->^BA) 

where 7^4 = 2 — 5/2(1 — g). They found that the 
case for the failure of RC was not proven within 
the limits of accuracy of the simulation, contrary 
to some earlier claims. 

They also considered Ising matter with a dl/l 
measure and found that the exponents were def- 
initely Onsager. It was found, for instance, from 
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the finite size scaling of the susceptibility maxima 
that 

- = 1.745(6) 

V 

so we can say with some certainty that 2D RC 
fails to reproduce gravitational effects when mat- 
ter is coupled in. Whether this is simply a failure 
of the scheme used in the coupling rather than 
the RC itself remains to be seen. 

The discussion of the RC results leads on to the 
broader question of what constitutes the univer- 
sality class of 2D gravity. There have been various 
results for Ising and Potts model simulations that 
have a bearing on this: 

• | |3^ ] Ising spins living in a flat geometry 
but with fluctuating connectivity, in effect 
a mixture of DT and RC approaches ^, give 
KPZ/DDK exponents. 

• js^] Ising spins on quenched ensembles of 
Poisonnian random lattices give Onsager 
exponents. 

• | |3^ ] Spins on quenched ensembles of 
2D gravity graphs may give "quenched" 
KPZ/DDK exponents. 

• js^] Squeezing 2D gravity (truncating the 
neighbour distribution in a triangulation to 
have just 5, 6, 7 neighbours) does not affect 
the KPZ/DDK values of the exponents. 

• | |3^ ] Spins on any, not just planar, cj)^ graphs 
(the duals of triangulations) give mean field 
exponents. 

In the above list appears to show that 

it is fluctuating connectivity that is important 
rather than curvature fluctuations in obtaining 
the KPZ/DDK exponents, whereas ||3^ suggests 
that even with a quenched ensemble of graphs, 
and hence no fluctuating connectivities, gravita- 
tional effects are still present. On the other hand, 
any collection of random graphs does not guar- 
antee gravitational effects: Poissonian random 



^The idea of a mixed DT/RC approach has also been sug- 
gested by Shamir in the 4D context and investigated ana- 



graphs give flat space jS^] exponents and generic 
(j)^ graphs give mean field exponents Q . 

Nonetheless, the vertex number distribution of 
2D DT triangulations can be quite brutally trun- 
cated and the effects of gravity still observed . 
Finally, the Onsager exponents seen in Regge cal- 
culus suggest that the fractal structure of the 2D 
DT lattices (apparently absent in RC, but still 
present in | |3^ , p7[ | ) may be an important factor in 
determining the universality class. Exactly what 
determines the universality class of the 2D gravity 
exponents is thus still not clear. 

We close this section with a quick word on 
plaquette surface simulations. Such models are 
closer to original Nambu-Goto idea in which the 
target space is discretized 



j dV-y/det \daXdbX\ 

so strictly speaking gravity has again disappeared 
from the picture. For a 2D surface embed- 
ded in 3D the general action S = (is Area + 
Pi Inter sections + (3cbends maye be written as an 
Ising like model 



ising 



lytically [^|. 



-I-J3 ^ aiGjUkOi 

where the Ising couplings are related directly to 
the surface couplings Ji = -|- /3c, J2 = 

"f-%, h = -f + x- Following this ap- 
proach for the Savvidy gonihedric string action of 
jl^ ] gives an Ising model with interesting proper- 
tics, including an unusual semi-global symmetry 
and Onsager- like exponents |^ . It has been sug- 
gested very recently Q that the gonihedric ideas 
that went into this string action may also be of 
relevance for 3D and 4D DT models, which we 
now move on to discuss. 

4. HIGHER D 

Some features of discretized 3D and 4D grav- 
ity in the DT approach appear rather similar to 
2D. There is a baby universe fractal structure and 
it seems that various scaling distributions can be 
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defined |^ in both the crumpled and branched 
polymer like phases. One feature of the 4D model 
not shared by its 2D counterpart was pointed out 
in p2| : The vertex order distribution in the crum- 
pled phase is very strange. The distribution p(o) 
of the vertex orders o{v), display a rather singular 
behaviour. 



K_2 = 2.0 
K 2 = 0.0 



P(o) 




o(v) 



10000 



In the branched polymer phase at K2 — 2.0 
there is a smooth distribution, as one might have 
naively expected, but the crumpled phase at K2 — 
0.0 has a detached "spike" at very large vertex or- 
der. Extensive simulations showed that this was 
not a finite size effect, it was not a thermaliza- 
tion problem and that it also appeared for other 
topologies. 

These observations were elaborated on in ||4^ 
where it was demonstrated numerically that the 
generic triangulation of the D sphere for D > 3 
contains one singular D — 3 simplex (k2 — for 
simplicity, but this is true throughout the crum- 
pled phase). This singular structure is attached 
to an extensive fraction of the rest of the trian- 
gulation. The geometric picture is shown above 
right for 4D and 5D respectively. 

In general one finds that the local volume asso- 
ciated with the singular D — 3 simplex grows as 
~ -^i^gj-eas tjj^g local volumes associated with 
the secondary simplices attached to this grow as 
~ V. Some counting arguments were presented in 
as to why this should be so. It is surprising, 
and perhaps rather disturbing, that such singu- 
lar features are present in the theory, especially 




given that the observed phase transition in 4D 
is apparently tied up with their appearance and 
disappearance. 

The nature of the 4D transition in DT has also 
been the focus of further work using the fractal 
blocking MCRG scheme. As in 2D it is the dis- 
tribution of geodesic distances that is of interest, 
which in discretized form is 



< r >= 



N2 

ij I 



In 4D, setting K2 = k and A^4 = N for concise- 
ness, one has to take account of the running of 
the coupling under the blocking move 

Sr — ri^SlnN + Vi^Sk 

so with the physical volume Na^ fixed, with a the 
lattice spacing, we have ^ In ^ = jSlnN. We can 
thus deduce Sn from measurements and extract 
the P function 



Sk 
(51ni 



The picture that emerged from doing this was of 
an ultraviolet fixed point [|6| , not an infrared one 
as had been suggested by some continuum treat- 
ments. One puzzling feature of the simulations in 
p6[ | was that the slope of the f3 function at the 
fixed point did not appear to be consistent with 
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a second order transition, which was the gener- 
ally accepted picture for DT in 4D. Further work 
looking at the scaling of cumulants such as 




where A'o = — + 2 was carried out to check 
this [Q. The finite size scaling form for this is 

C2(iV4)^A^4/((«2-'*2)^4) 

where b — c ~ 1 for a first order transition and 
take non-trivial values for a second order transi- 
tion. The simulations found the first order values 
b = c ^ 1. Further evidence for the first order 
nature of the transition was found by deBakker 
]45| who saw time series flips characteristic of first 
order transitions, as well as a double peak his- 
togram structure emerging on larger lattices. 

Simulations discussed at this conference Q 
show that the first order transition is not an ar- 
tifact stemming from the inclusion of the "/{N^ — 
y)^ term in the action - relaxing the constraint 
to allow huge volume excursions has no effect, so 
we are left with the task of deciding whether the 
model should be patched up to give a second or- 
der transition, is still useful as it stands, or should 
be consigned to the bin. 

The optimists argue that some things are right 
and in any case there are lattice theories with 
first order transitions and a Coulomb phase for 
an entire range of couplings. The power law cor- 
relations in the branched polymer phase might 
be argued to be just such a case |Q 0. The DT 
model also appears to display gravitational bind- 
ing ^ as defining one and two particle connected 
propagators for scalar test particles shows 

If we insist on getting a continuous transition, 
various alternatives suggest themselves. Ray 
Renken presented evidence using a general- 
ization of the node decimation scheme to higher 
dimensions that suppressing vertices of high or- 
der with a measure term could weaken the 3D DT 
transition. A similar effect may also appear in 4D. 

tWe should remember the caveat of the polymer example 
|13| , however. 

*4D Regge Calculus would also appear to be attractive 

i- 



It might also be possible to use the ideas of |40| 
stemming from the gonihedric string, where new 
integral invariants for simplicial complexes corre- 
sponding to higher derivative terms were derived. 
These would have the effect of "stiffening" the 
lattices and possibly changing the nature of the 
transition. One might, of course, lose the transi- 
tion completely with either of these approaches. 

If we are willing to overlook the problems expe- 
rienced in 2D by Regge Calculus as an aberration 
of lower dimensional gravity alone it is interest- 
ing to ask what what Regge calculus says about 
the 4D transition. This is not such a leap of faith 
as it might appear at first sight - in 2D there is 
no classical action at all and it is the quantum 
effects that contain the dynamics of the theory, 
whereas in 4D there is a non-trivial classical dy- 
namics. It is thus conceivable that a model like 
RC which starts out "close" to smooth manifolds 
might do a better job in 4D. The Vienna group 
have looked at various models in 4D including 
an Ising-link model in which the edge lengths are 
restricted to two values qi = bi{l + eai) where e 
is a small constant and ai is an ising spin living 
on the link. They found a first order transition at 
positive /3 (Newton coupling) as well as another 
second order transition at negative (3. 
If one couples in SU{2) gauge matter 

S = Atat - ^ ^ W^t Re[Tr(l - {/,)] 

t t 

the gauge field string tension scaling at the neg- 
ative (3 transition is still what would be expected 
of a physical theory. The suggestion is then that 
somehow DT is looking at the "wrong" first order 
transition - the physical theory being the one that 
resides at the negative (3 second order transition. 

5. IN CLOSING 

In summary: the ID results shed an interest- 
ing light on polymer phases in higher dimensional 
models, 2D DT is in very good shape, 2D RC less 
so, in 4D given the recent results on the order of 
the transition the big question is whether either 
the DT or RC theories are telling us anything 
about Euclidean quantum gravity. There will no 
doubt be much of interest to report at Lattice97. 
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Lest the reader be too despondent about the 
implications of the 4D results I would like to fin- 
ish with a couple of sentences on recent purely 
analytical work in Ashtekhar's approach ||5^ to 
quantum gravity that suggest a model not far re- 
moved from Regge Calculus or DT emerging from 
a completely different perspective. If we make 
3-1-1 split, and use the canonical variables 



3. 
4. 



EfE" 



where the 3-metric, Ef a dreibein and K 

the extrinsic curvature, then a canonical transfor- 
mation takes us to Ashkekhar's "new" variables 

where G is the Newton coupling. The constraints, 
a major stumbling block in the canonical ap- 
proach for many years, are greatly simplified with 
these variables at the expense of a complex con- 
nection A. 

As in lattice gauge theory a loop representation 
exists and it was recently realized that a spin net- 
work basis resolves the Mandelstam constraints 
that had frustrated progress in this approach. 
The picture that one arrives at is, crudely, of a 
embedded trivalent |^ graph with edges coloured 
by spins j such that |ji - j2| < js < Ji + J2- 
We can even have real connection A if we can 
put up with more complicated constraints (i.e. 
Hamiltonian) ||5^ . Continuum spacetime emerges 
from coarse graining this lattice-like structure. To 
quote Ashtekar js^]: "The fundamental Planck 
scale excitations of the gravitational field arc 1- 
dimensional and the corresponding geometry is 
distributional" , so perhaps the next presenter of 
this review will not be forced to include "on the 
lattice" in the title if space(time) is a lattice. 
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